Effective Hamiltonian study of excitations in a boson- fermion mixture with 

attraction between components 
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An effective Ifamiltonian for the Bose subsystem in the mixture of ultracold atomic clouds of 
bosons and fermions with mutual attractive interaction is used for investigating collective excitation 
spectrum. The ground state and mode frequencies of the *^Rb and mixture are analyzed 
quantitatively at zero temperature. We find analytically solutions of the hydrodynamics equations 
in the Thomas- Fermi approximation. We discuss the relation between the onset of collapse and 
collective modes softening and the dependence of collective oscillations on scattering length and 
number of boson atoms. 
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I. INTRODUCTION 

Bose-Einstein condensation (EEC) in ultracold 
trapped atomic gases [ll-Q has been the subject of 
intense theoretical and experimental interest [1, 
Experimental studies of the BEC properties in confined 
vapours of alkali atoms have been extended to double 
Bose condensates [B-Q, to the achievement of quantal 
degeneracy in gases of fermionic atoms llOl. and to 
dilute mixtures of Bose and Fermi particles |lll - [l5l | . 

CoUisional interaction between bosons and fermions 
greatly affects the properties and leads to a rich phase di- 
agram of the degenerate mixtures. Theoretical consider- 
ations predict the phenomenon of component separation 
for systems with a positive coupling constant il6H2Q] , in- 
stabilities and significant modification of the properties 
of individual component in a case of boson-fermion at- 
traction [2l| - |23j . the formation of a superfluid state due 
to boson- induced fermion- fermion attraction [24j] . The si- 
multaneous collapse of the two species has been observed 
experimentally in the ^"K - ®^Rb mixture by Modugno 
and co-workers [T^ as a sudden disappearance of fermion 
cloud when the number of bosons is increased over an in- 
stability value Nbc ~ lO'^, the critical number of fermions 
being TV^ ~ 2 x 10". The stability diagram for the 
- ^^Rb mixture and the critical particle number for the 
onset of the collapse has been considered in [25j . 

The dynamical properties have been investigated for 
various boson-fermion mixtures, i.e., '^^K - ""^K ^6, 27], 
^Li - ^Li dldl, 87Rb - [33. The theoretical studies 
of the collective oscillations based on a random-phase ap- 
proximation [2^, [2^ [2^ , direct numerical integration of 
the time-dependent equations [13] , semi-analytical meth- 
ods [m, [13, or on numerical solutions of the coupled 
eigenvalue equations for boson and fermion density fluc- 
tuations [305. 

The analytic treatment of the elementary excitations 
for Bose- condensed gases confined in magnetic traps has 
been given by Stringari [ss!]. The Stringari's dispersion 
law of the discretized normal modes does not depend 
on the interparticle interaction strength and holds in 
the hydrodynamic limit Nbub / o^ho ^ 1, where a^o — 



{h/niBi^By^^ is the harmonic oscillator length, qb > 
is the bosonic scattering length, and Nb is the number 
of bosons. The numerical confirmation of these behavior 
has been given in [s^l by applying linear-response theory. 
In the case of attractive interaction (as < 0) the system 
becomes more compressible when approaching the criti- 
cal number Nbc for collapse. The Nb dependence of the 
lowest monopole frequency uj^ can be determined ana- 
lytically [35] using a variational calculation of the ground 
state based on Gaussian trial wave function for the order 
parameter. It vanishes as ujm ~ (1 — Nb/Nbc)^^^ ■ 

Collective modes of a boson-fermion mixture with mu- 
tual attractive interaction have been studied by Capuzzi 
and coworkers [sS] using the equations of generalized 
hydrodynamics. The numerical procedure included the 
decomposition of the density fluctuations into compo- 
nents of definite angular momentum Z, i.e. 5p{v) — 
5pi{r)Yim(0,if), and solving an eigenvalue problem for 
coupled equations in a given I subspace by means of stan- 
dard linear-algebra routines. It was found that the col- 
lective spectra show a frequency softening of a family of 
modes of bosonic nature as a signature of the incipient 
collapse. This softening becomes most pronounced in a 
very narrow region of parameters near collapse. 

The purpose of the present work is to present an an- 
alytical method for studying the dynamics of a harmon- 
ically confined Bose- Fermi clouds and provide a theo- 
retical discussion of interaction effects on the collective 
oscillations. We analyze quantitatively dynamical prop- 
erties of the ^''Rb and ^^K mixture with an attractive 
interaction between bosons and fermions at T = . The 
dynamics of a boson-fermion mixture is found amenable 
to analytic approaches by effective Hamiltonian method 
[sgI . [13] ■ This method well describe the density pro- 
files of the Bose- Fermi mixtures even close to collapse 
and predicts the correct value for critical particle num- 
ber to make the collapse diagram [25]. In order to gain 
more physical insight and to solve analytically the lin- 
earized hydrodynamic equations we consider a special 
limiting case in which the kinetic energy term is negli- 
gible compared to confining and boson- boson interac- 
tion energy, the Thomas- Fermi (TF) approximation for 
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equilibrium condensate profile similar to Stringari's ap- 
proach [sl]. But in contrast to [s^ our spectrum explic- 
itly depends on the ratio of the boson- boson and boson- 
fermion couplings. We show a signature of the softening 
of the bosonic collective mode spectra in the vicinity of a 
collapse, resulting from peculiar property of the ground 
state density when the system become more and more 
compressible in the vicinity of the collapse. 

The paper is organized as follows. In Sec. II, we 
present our effective Hamiltonian approach and intro- 
duce necessary notations. In Sec. Ill we give our results 
for ground state boson density distribution and collective 
mode spectrum. We then summarize in Sec. IV. 



II. THEORETICAL MODEL 

A. Effective interaction 

Our starting point is the functional-integral represen- 
tation of the grand-canonical partition function of the 
Bose- Fermi mixture. It has the form [s^, |3^ : 

Z = J D[r]D[<j>]D[r]D[^j]expl^-^[SB{<l>*,c^)+ 

+ SF{r,^j) + s,„t{r,<i^,r,m, (i) 

and consists of an integration over a complex field ^(t, r) , 
which is periodic on the imaginary-time interval [0, 
and over the Grassmann field iPItjT), which is antiperi- 
odic on this interval. The field (j){T, r) describes the Bose 
component of the mixture, whereas iP{t, r) corresponds 
to the Fermi component. The term describing the Bose 
gas has the form: 



Sb{4>\4>) = j^^ dT J dr|0*(T,r) (^h^ 



2m B 



+ 



Because s- wave collisions between fermionic atoms in the 
same hyperfine state are forbidden by the Pauli principle 
the Fermi-gas term can be written in the form; 



dr 2m p 
+ V^i.(r)-Mf)V(r,r)}. (3) 



Here VB^piy) are the external confining potentials, and 
IJ.B,F are the chemical potentials for the Bose- and Fermi- 
components respectively. Under isotropic harmonic con- 
finement Vb,f(j) = "T-B.Fi^B f'^^/^i '^b and mp are 
the masses of bosonic and fermionic atoms respectively. 
The trap parameters ujb and ujp are chosen in such 
a way that mB^Js/'^ ~ mpLUp/2 — Vq, that is why 
ujp = yJmB IrnpUJB ■ The radius of the Fermi- cloud R p 
can be estimated as VqE?p ~ /if- 



The term describing the interaction between the two 
components of the Fermi- Bose mixture is: 

S,„t{rA,r,i')^9BF J drj dr|^(T,r)n</,(r,r)|2, 

(4) 

where qb — 'iirh'^aB /mB and gBF — 2Trh^aBF /mi, 
TO/ = mBmp I [mB -t- mp), gb and ubf are the s- wave 
scattering lengths of boson-boson and boson-ferniion in- 
teractions. 

The integral over Fermi fields 



Zp 



D[i!*]D[ilj]e 



(5) 



is Gaussian, we can calculate this integral and obtain the 
partition function of the Fermi system as a functional of 
Bose field (^(r, r) which for T = has the form [H, H^]: 



Zr. 



cxp 



hp 



Seff). S^ff^ / dr / dr/e//(|0(r,r)|) 



/e//(|0(T,r)|) ^--y,[y,p- Vp{v) - aBp\4>{T,v)\^] 



5/2 

(6) 

where x — \/2m'l^^ / (in'^hP) . 

Using the fact that due to the Pauli principle (quantum 
pressure) the radius of the Bose condensate is much less 
than the radius of the Fermi cloud i?F ~ \/ Hf/Vo, one 
can use an expansions in powers of Vp(r)/ up and obtain 
the effective Hamiltonian in the form |36l . |37[ : 



Z = 
hp 

dr 



d^rcj)*{T,r)h 



d 







+ VB{v)-^,B)<t^{r,v) + ^-\ct^{r,v)t]. (2) ^^//I'^*.'^] 



dr 



dr 

re 

2tob 

2 



5b4 



Seff = 



(7) 



+ (t4//(r)-A^B)|0| 
where 



BE 
9eff 



(8) 



; ^^ 1/2 \ 

Ko = (1 - 2^A^F 9BF), 



^BB 



1/2 „2 



9eff = 9B - 2^/^F 9bf 



9eff 



3>€ 



9%F- 



(9) 



The first three terms in ([5]) have the conventional Gross- 
Pitaevskii [ioj form, and the last term is a result of boson- 
fermion interaction. It corresponds to the three-particle 
elastic collisions induced by the boson-fermion interac- 
tion. In contrast with inelastic 3-body collisions which 
result in the recombination and removing particles from 
the system [41], this term for qbf < leads to increase 
of the gas density in the center of the trap in order to 
lower the total energy. 



3 



B. Hydrodynamic approach 

Now rewrite the action in terms of hydrodynamic vari- 
ables density and phase (j){T,r) = ^yp{^yr)e^^^^''''' 



Z = J D[p]D[0] e-*(^«+^=//), 



Sb 



hi3 

Seff ^ I dr I d 



Srrif 



/'V7 \2 r,BB „BF 



(10) 



To simphfy the formahsm we introduce dimensionless 
variables rescaled by the natural quantum harmonic os- 
cillator units of length a ho = \/ fi/'m-Bi^B, and energy 
hujB- r = ahoTc', T = t' /uB, E = HujbE', p p' la\o^ 
P = hujBP'/al^, S = S'/h. Then the effective Hamiho- 
nian takes the form (the primes omitted) 



Taking t he g radient we obtain quantum Bernoulli's equa- 
tion l40l.l4al 



9v 

at 



= 0. 



(17) 



Equations ()13p . ()17p are the modified hydrodynamic 
equations of Bose- fluid, which incorporates the effect of 
Bbose- Fermi interaction. They correspond to the equa- 
tion of state of the mixture where pressure and density 
are related by 



P = -up' 



^vp- 



(18) 



where Pkm (p) is a kinetic energy contribution associated 
with spatial variations of the condensate density p{r). 

Now to obtain elementary excitations one linearizes 
Eqs. ([T^. PT|) around ground state solution /9o(r), vo = 
and looks for eigenvalue mode Sp{r,t) = Sp{r)e~^'^*. 
The hydrodynamic amplitudes can be combined into a 
sing le second- order equation for the density perturbation 



Seff = I dT 



Heff[p,e] 



,s 90 

^ '^^P-Q:;: + Heff[p, 



IjVpf 
8 p 

2 , «„3 



+ (K//(r)-AiB)p+^p'+ gP' 



where K//(r) — A;or^/2, and we introduced dimension- 
less parameters for boson- boson and boson- fermion cou- 
plings u = 2g^/j/al^fujB and v ^ 2gf/^/4„fiwB. 

To obtain equation of motion of Bose condensate den- 
sity and phase one varies the action Sb + Seff and goes 
to the real time t — it [s^. Variation over the phase 



SiSE 



0, 



dt 59 



0, 



yields the continuity equation 



|^H-div(pv)=0, 



(12) 



(13) 



where v = ft/TO[V0(r, t)]. Variation over the density 
p(t, r) 



5{Sb + Seff) 



= 0, 



SH, 



Sp ' dt Sp 

yields the equation of motion for the phase 



= 0, 



d9 v2 
dt^Y 

72 



„2 



^ = - ^^ + '^°y + 2^+2^ - 



uj'^Sp + div(poV(5/i) = 0. 



(19) 



Solutions of give the low- frequency condensate col- 
lective modes of an inhomogeneous Bose component with 
a local condensate density /9o(r). The Eq. (ITOl) replaces 
two coupled eigenvalue equations for the density fluctu- 
(11) ations SpB,F of each species [H, [s^l- The coupled eigen- 
value equations predicts two sets of eigenvectors [28l.l30j. 
which can be labeled as fermionic and bosonic ones, ac- 
cording to the nature of their eigenvalue in the limit of 
vanishing gBF- Solutions of (IT9t describe oscillations of 
the bosonic cloud on the background of the fermionic 
component and corresponds to that branch of excitations 
which has the bosonic origin. 



III. RESULTS AND DISCUSSION 

A. Ground state 

To clarify the main features of excitation spectrum of 
a Bose- Fermi mixture we illustrate our results on the 
example of *^Rb - '^''K mixture and consider an isotropic 
trap when the problem can be treated effectively as a one- 
dimensional. In the experiment of Modugno et al. (T^ 
K and Rb atoms experience potentials with an elongated 

(14) symmetry with substantial value of trap asymmetry pa- 
rameter A = ujz/uJ_L. To compare our results with those 
in the experiment we have to rescale number of bosons 
Nb by the reverse trap asymmetry ratio Nb — > Nb/X- 
For the chemical potential of an ideal Fermi gas in a trap 

(15) one can use the relation fj,p = fujjFiQXNpY^^ j43l |. 

The parameters of the ^^Rb and mixture are the 

(16) following [ll|: Np = 2 ■ lO'', Ub = 5.25 nm, ubf = 
— 21.7^4 g nm. The magnetic potential had an elongated 
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FIG. 1: The profile of the ground state condensate wave func- 
tion at ric — 20 and nj. = 10 (solid line). The line with circle 
markers corresponds to the ideal Bose gas in an external har- 
monic confinement. Dashed line is for the density profile in 
the TF approximation (|22p . 



symmetry, with harmonic oscillation frequencies for Rb 
atoms uj± = lob ~ 2tt x 215 Hz and u}b,z = Awb = 
27T X 16.3 Hz. At this parameter values the reverse trap 
asymmetry ratio 1/A = 13.2, characteristic length a^o — 
735 nm, the chemical potential for fermions fip ~ 31 ftws, 
ujF « lAToJB, ko = 1.07, u = 0.11, v = -0.0003. 

The ground state density distribution po{r) is defined 
by the stationary equation v = 0, Jl = /is and gives rise 
to the modified Gross- Pitaevskii equation 



2^ 



(20) 



The condensate density p(r, t) is normalized to the 
number of atoms in the condensate J (Prp{r, t) = N . 
In the T — > limit, N coincides with the total number 
of bosonic atoms in the trap. 

Ground state properties of the mixture and numerical 
solutions of Eq. PU)) have been considered in [3, l45| . 
It was shown that the TF approximation gives the good 
description of density distribution and BF mixture can 
be accurately considered in the TF approximation (46| . 
In TF approximation 



k^r"^ - 2fiB + upa + vpl ^ 0. 



In this case the density profile has the form 



Rr 



Rc 



r ^ Rb, 



(21) 



(22) 



where Ucr = u/{2\v\), R^^ = u^/(4|?;|), Rb = R/Vk), 
i?^ = 2^B, = Rl^ — R"^. In the limiting case of non- 
interacting Bose and Fermi clouds {gBF 0, v 0, 
a ^ 1, Rcr ^ a + R^ /2a, ko « 1, ncr/2Rcr — we 
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FIG. 2: Chemical potential pB (in the unit of /kjs) and num- 
ber of bosons Nb as functions of the rescaled central density 
n'r- 



recover the TF distribution for single bose- condensate 
Po(r) = (i?2 - T^)/u. 

The parameter i? up to a multiplicative factor 1 / ^/ko 
is the radius Rb of the Bose condensate. In TF ap- 
proximation Eq. ((2T|) for the center of the trap enables 
to relate the value of R, i.e. /i, with central density as 
R^ = upo + vpq. Let us express R through parameter n'^ 
as 



R^ = <(1 - 



(23) 



The values of the parameter R gains the maximal value 
R = Rcr at < = u^/2\v\ w 20. 

The numerical solution of Eq. ((20|) for the condensate 
wave function 0(r) = ^/poir) and TF profile ([^ are 
presented in Fig. [T] They parameterized by rescaled 
central density = uric, Uc = |0(O)p. Note that the 
wave functions in the figure are normalized to unity. 

Figure [T] shows two characteristic profile of the con- 
densate wave function with increasing the central den- 
sity from n'^ = 10 to n'^ = 20. For n'^ < 20, the behavior 
is characteristic of the Bose gas with repulsion, i.e. the 
evolution of the profile corresponds to a monotonic ex- 
pansion of the boson cloud with increasing number of 
bosons. The cloud density becomes more flat at the trap 
center, approaching TF analytical solution. For n'^ > 20, 
the solution changes qualitatively: the central density 
begins to increase. 

For n'c > 20, the central density increases significantly 
in a small region near the trap center. We attribute the 
solution with n'^ > 20 to the nonstationary states of the 
condensate through which the collapse of the condensate 
wave function occurs. The value « 23 corresponds the 
critical particle number Nbc — 8 • 10^ above which col- 
lapse occurs. In that point the branch of stable solutions 
meets the branch of unstable solutions of the GP equa- 
tion. This qualitative behavior is the generic signature 
of a Hamiltonian saddle node bifurcation [47j . 

Fig. [2] shows the boson chemical potential pB and the 
number of bosons Nb as function of the parameter n'^, 
obtained at numerical solution of Eq. d^H]) . Both curve 
have maximal value at n'c — 23. The condition dN/dn'c = 
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is a sufficient condition of the zero excitation mode 48^ . 
It means that it does not cost energy to deform the initial 
density profile continuously into the final, which indicates 
the threshold for instability 



d^rpoir,n'J = N{n'^). (24) y{l-y)w"+(^ 



radial eigenfunction 5pni{r) and a spherical harmonic 
YimiO^ip), where n is the radial quantum number. The 
associated eigenfunctions are polynomials of the form 

Spniir) - Cr'«;„i(rVi?2) Q{R ^ r), (29) 

where Q{R — r) is the step function. The polynomials 
Wni (y) are polynomials of order n and satisfy equation 

2/ + 3 21 + 5 \ , e^l 

—w = 0, (30) 



y u) 



It is interesting to note that the radius of the con- 
densate Rb = R/y/ko for TF function (f23|) qualitatively 
conform these behavior, i.e. it increases for n'^ < n^g and 
decreases for n'^ > ?ij,g. There is a critical point where 
the dependence Rb{pb) shows the maximum. This re- 
markable feature indicates the onset of a collapse in the 
system [3]. So we can anticipate some features in the 
dynamics of the bose- cloud near the collapse and in the 
eigenfrequency spectrum. 



B. Collective excitation spectrum 

We now give a brief overview of the calculation of the 
low- frequency collective mode spectrum. We have to 
give explicit solution of Eq. (fT9|) . The main difficulty 
associated with Eq. (fT9|) is the proper treatment of the 
differential operator term div(poVi5/I) with 

<5/I--^5(^-^ AVp^ +^5p + i;p(5po. (25) 

For uniform Bosc fluid with constant condensate density 
Po the first term accounting for quantum pressure re- 
duces to — l/(4po)A5p and the differential operator term 
is handled in the straightforward manner 



div(poV(5p) Po A(5/i 



A?5p 



upa 



ASp. (26) 



Eigenvalue modes are the plane waves Spgjr) = Ce^'^^ 
and the resulting eigenvalues are given by [51] Wq = 

wpoqV2 + qV4. 

For an inhoniogeneous Bose- condensate in a trap the 
TF approximation must be called on to describe ana- 
lytically the condensate collective modes [s^. For the 
ground state the TF density has the form 



i?2 



r^R, 



(27) 



where the radius of Bose- condensate is i? = 2p,B- In Eq. 
(1251) one omits the quantum pressure term, 5]! ~ u5p/2, 
and 



where e = w^. They can be expressed through the hyper- 
geometric function Wni{y) = F{a, f3,^,y) with parame- 
ters a — n, (3 = I + 3/2 + n, = I + 3/2. The associated 
energy eigenvalues are found to be 



uli ^ ujl{2n^ + I + n{2l + '3)). 



(31) 



The hypergeometric function F{a, I3,^,y) satisfies the 
equation [52| 

y{l - y)F" + [7 - (a -f /3 + l)y]F' - a(3F = 0. (32) 

The eigenfunctions for the Bose-condensed cloud 

fluctuations vanish outside the cloud radius R and 
present a discontinuity at R. The discontinuity is phys- 
ically acceptable in view of the fact that the kinetic en- 
ergy term has been set as negligible in taking the strong- 
coupling limit. 

Let us now address the issue in the case Bose- Fermi 
mixture. As well as in the previous case we omit the 
quantum pressure term in (j25p and the perturbation for 
Jl takes the form 



f^M = -j^Sp + vpoSp, 



and 



^5p=[-+ wpoj V(5p + t;(Vpo)(5p, 

(33) 



div(poV(5p) = Qpo + f Po) A(5pH 



^Vpo + ^wVpg ) -ySp+vpl 



\ Po 



Apo 

Po 



5p. 

(34) 



Due to Vpo/po ^ 1 we should omit the last two terms in 
square brackets in TF approximation. So we can rewrite 
the expression (j34p in the form 



div(poV(5p) 



1/ 2n 2 

-(wpo + i;po) + -Po 



ASp+ 



\ {upo 



vPo) 



VSp. (35) 



Now we use in Eq. (1351) the ground state density profile 
Po (r) from Eq. (|22|) , which simply states that due to Eq. 
(PT|) upo + fpo = 2pB ^ k^r"^ and come to the result 



div(poV(5p) 



-A5p 



R' - r ' 



VSp (28) 



For a spherical trap, an eigenfunction of Eq. ([T^ can 
be written as a product Sp{r) = 5 pni{r)Yira{0 , of the 



div(poV5p) 



i?2 



i?2 



2 

fcor^ 



A5p' 



vpl 



VSp. (36) 
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Note that in the — > limit the Eq. (|36|) transforms 
into Eq. (|28p for single Bose condensate. In principal, 
the solution of Eq. (|19p with approximate div(poV(5/i) 
giving by Eq. p6p should give rise to collective excita- 
tions spectra in TF approximation. But the term wpg is 
the complicated function with square root singularity. In 
order give analytical treatment of collective excitations 
we have to approximate the wpp term by more analyti- 
cally simple function. It is reasonable to simplify the vp^ 
term and approximate it by quadratic function on the 
radius r. For we use Eq. ((22|) and expand the radical 
near the center of the trap which gives rise to 



1.2 



(37) 



where R'g = ab/ko and b = Rcr — a. In the framework 
of the density profile (|37|) it is possible to give explicitly 
solutions of Eq. (|19p through the hypergeometric func- 
tions. 

It is clear that the result ([57]) does not describe den- 
sity near r = R'^. This fact is an artifact of the 
TF approximation, arising from the vanishing of po(r) 
at r = R'g. The dependence Pq in Eq. (|37l) has 
slightly different radius of the cloud R'g in comparing 
with Rb in Eq. (j22p . but this approximation is consis- 
tent with the TF approximation because in the region 
near boundary density of the cloud is small. Note that 
(i?^)2 ^ R-i/ko - (i?2^ - ^i?2^ - R^)/ko and R'g ~ Rb 
at i? ^ Rcr- In this regime the density profiles giving 
by Eqs. (1^^ and ([57]) are the similar. In contrast in the 
critical regime at i? < Rcr the parameter R'^ <^ Rb and 
the approximation Eq. ([37)) is no longer valid. 

Now for vpq term in Eq. (p6| we use expression (|37|) 
and introduce notations ai = ko{l — b/a), oli — fco(l — 
2b I a). Then 



div(poVp^) = ^^^^^^ASp - a,ri-Sp. (38) 
Now the equation takes the form 



d_ 



i?2 - d 

^L^Adp~a3r—Sp + e5p = 0. (39) 
Z or 

where R^ — k^R^g/ai, a^, — Q2/«1i e — uP' jcti- Separat- 
ing variables 8p(v) — r'' G{r)Yim{d,(p) and substituting 
y = r^/Rl^ one obtains for the function w{y) — G{r{y)) 
the following equation: 



21 + ?, 21 + 2,- 



y{l-y)w"+[ ^''~^'' - ""^^^ "^'^ yj w'h = 0. 

(40) 

If we compare it with analogous Eq. (|30|) we can see 
that it contains additional parameters ai and a2 includ- 
ing the effect of interparticle interaction. At u — > the 
coefficients ai,a2 — >■ 1, and the Eq. (|40l) transforms into 
Eq. ([50)) for a pure boson system. Eq. ([HH) is a standard 
equation of the form ()32[) for hypergeometrical function 
F(a, P^"f,y). It should be a polynomial of n- th degree. 
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FIG. 3: The coefRcients ai and 02 of the formula (|4ip as func- 
tions of A'^B and |a_BF|. Left panel: u = 0.11, v = —0.0003. 
Right panel: u = 0.11, iVs = 5 ■ 10^. 



which for v — ?> tends to Stringari's solution ([^^ . So the 
parameters of the hypergeometrical function F{a, /3, 7, y) 
get the values a — —n, (3 — {21 + 1 + 2a-^l2 -f n, and 
7= (2/-h3)/2. 

The associated energy eigenvalue are found to be (in 
dimensional units) 

= wl [ai(2n2 + n(2Z + 1)) + a2{2n + V)\ . (41) 

The result ([4T)) can be considered as a qualitative sig- 
nature of shifts of eigenfrequencies in the presence of 
the boson- fermion coupling. At z) — t' frequencies ([IT)) 
transform into the Stringari's spectrum (I5T]) . 

The dispersion law of the normal modes given by the 
formula (HTI) has the additional coefficients a\ and a-i 
which incorporate the effect of interparticle interaction. 
In explicit form the coefficients a\ and a2 are related 
with Nbi Np, and qbf through the parameters R and 

Rcr 

( Rcr-VRlr^R' \ 

ai ~ ko \ I = 

V ^Rl-R^ J 

^^ Rcr-VRl-RA 



0.2 = k 



(42) 



Formally these relations are valid only at i? ^ Rcr- In 
our approach we interpolate they up to i? < Rcr- These 
rather crude approximation nevertheless gives rise to the 
qualitatively correct picture of collective excitation spec- 
tra tendency to become softening near the collapse tran- 
sition. 

The coefRcients a\ and are plotted in Fig. |3|as func- 
tions of the number of bosons Nb and scattering length 
IobfI- The chemical potential /is, and the number of 
bosons Nb are fixed by Nb — J Po(r)c?'^7'. For small 
values of the IobfI and Nb the coefficients ai and a2 
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FIG. 4: Frequencies (in units ujb) of low- lying monopole 
modes Z = 0, n = 1,2,3 (solid line), dipole and quadrupole 
modes n — 0, I — 1,2 (dashed-dotted line) as functions of the 
number of boson atoms Nb- 
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FIG. 5: Frequencies lu„i (in units us) of low-lying monopole 
modes Z = 0, n = 1,2,3 (solid line), dipole and quadrupole 
modes {n — 0, I — 1,2) (dashed-dotted line) as functions of 
the scattering length ubf for Nb = 5 ■ 10"^. 



approach to one, while for Nb and IubfI large enough, 
one observes, as expected, important deviations due to 
the boson- fermion interaction effects. The figure shows 
that the system will be in the collapse regime both at 
sufficiently large Iubf] and Nb- 

We will further limit the discussion to the case of exci- 
tations of low multipolarity corresponding to the energy 
range ujni ^ where prediction (j4ip is expected to be 
accurate. To illustrate the dependence (HIT) we consider 
the low-lying monopole modes {I = 0) and the surface ex- 



citation modes (n = 0), and keep the number of fermions 
fixed. 

In Fig. |3]wc show the evolution of utni as a function of 
number of bosons for aBF = —21.7 nm. Three low- ly- 
ing monopole modes corresponding n = 1,2,3 and dipole 
and quadrupole modes are shown. At increasing values 
of the boson number the frequencies of modes display a 
monotonic decrease and ultimately vanish in the vicinity 
of the collapse. Points of a collapse are slightly different 
for different modes. This is a result of a rather crude 
approximation that we chose for the equilibrium density 
po(r)- The second source of errors comes from the TF 
approximation. Results based on the TF approximation 
are in fact not adequate in the vicinity of r ~ Rb and 
beyond. If one works with the full hydrodynamic the- 
ory keeping the kinetic energy contributions in (1191) , one 
finds that po{r) does not abruptly vanish at r = Rb but 
exhibits a small tail which slowly goes to zero. This tail is 
neglected in our approximation. Nevertheless the behav- 
ior correctly display the situation which occurs near the 
collapse. The similar frequency behavior of the bosonic 
modes as a signature of the incipient collapse has been 
obtained Capuzzi et al. with numerical calculations on 
the basis of generalized hydrodynamics equations [30|. 

Fig. [5] shows the frequencies w of the same low- lying 
monopole and surface modes as in Fig. 2] as functions 
of the scattering length ubf at fixed number of bosons 
Nb = 5 • lO'^. At larger values of jasFl the eigenfrequen- 
cies show softening as collapse is approached. In contrast 
with the previous picture the points of the collapse are 
approximately the same for different modes. At increas- 
ingly large boson- fermion attraction the densities of the 
two species increase in their overlap region and collapse 
occurs when this attraction overcomes the Fermi kinetic 
pressure and the boson- boson repulsion. Similar behav- 
ior has been revealed in ''Li-^Li mixture with numerical 
calculation of eigenfrequencies (30| . 

Note that both components of the *^Rb-''°K mixture 
are inside the same magnetic trap, but their trapping fre- 
quencies differ considerably as a consequence of the large 
difference in atomic masses. At this conditions the gen- 
eralized Kohn theorem [i^ is satisfied when the species 
are uncoupled {ubf 0) and each component enable to 
oscillate with its own frequency. From Fig. [S] one can 
see that the dipole mode frequency w/j (n = 0, Z = 1) in 
broad region on aBF except near the collapse regime is 
found to be not affected by the interactions and equals to 
the confining potential frequency ljb- Iu fact this mode 
corresponds to the oscillation of the center of mass of 
the Bose- gas with ujd = lob driven by the external har- 
monic potentials. The same one can say on the monopole 
ujm and the quadrupole ujq mode frequencies. In this 
limit (aBF ^ 0) one recovers the results ojm = V^i^B 
{n = 1,1 — 0) and ujq = V^ujb (n — 0,1 = 2) predicted 
by the hydrodynamic theory of superfluid [33i] . 

In Ref. [sil, Is^l the collective oscillations were calcu- 
lated by the sum-rule approach and coUisionless random- 
phase approximation (RPA). Mixing angle of bosonic and 
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fermionic multipole operators was introduced so that the 
mixing characters of the low- lying collective modes were 
studied as functions of the boson- fermion interaction 
strength. For an attractive boson- fermion interaction it 
was shown that the low-lying monopole mode becomes a 
coherent oscillation of bosons and fermions and shows 
a rapid decrease in the excitation energy towards the 
instability point of the ground state. In contrast, the 
quadrupole mode give indications of frequency increas- 
ing and it contradicts the behavior of modes in Ref. [13] ■ 
In Ref. [lOl instead it was found that, as the mixture 
is driven towards collapse instability, the frequencies of 
the modes of bosonic origin show a softening, which be- 
comes most pronounced in the very proximity of collapse. 
Explicit illustrations of these trends were given for the 
monopolar spectra. 

A comparison of the hydrodynamic spectra (30| with 
the spectra calculated on the basis of the effective Hamil- 
tonian suggests that eigenvalue equations with the dy- 
namical coupling between the two components yield the 
similar bosonic mode behavior both as functions of Nb 
and asF SiS the hydrodynamic- type equations (jl9p with 
the fermion degree of freedom being integrated out. Soft- 
ening of the bosonic mode spectrum is a general signature 
of the dynamics of the Bose- Fermi mixture with attrac- 
tion between components arising from peculiar properties 
of the ground state and instability point near the collapse 
transition. 



IV. SUMMARY 

We have studied the collective mode frequencies in the 
boson- fermion mixture. To this purpose we have used 



the effective Hamiltonian for the Bose system, where the 
fermion degrees of freedom are integrated out [H, [STj . 
The effective Hamiltonian incorporates the three- parti- 
cle elastic collisions induced by the boson- fermion inter- 
action. In terms of coupled eigenvalue density fluctua- 
tion equations for the bosonic and fermionic components 
it means that we exclude fermion fluctuations from the 
bosonic density fluctuation component resulting in more 
involved structure of bosonic fluctuations. 

This approach enables to account for interparticle in- 
teraction strength in collective excitation frequencies and 
can be considered as an extension of Stringari's solution 
[3^ for Bose - Fermi mixtures. The analysis of the mode 
frequencies as function of the number of boson atoms and 
the boson- fermion coupling strength yields a dynami- 
cal condition for a collapse of the system. The behavior 
of mode frequencies is consistent with numerical analy- 
sis on the basis of generalized hydrodynamics equations 
[30| . The effective Hamiltonian approach proves to be ef- 
ficient tool for the analysis of both the ground state and 
dynamical properties of boson- fermion mixtures and can 
be considered as alternatives to Capuzzi's et al. [13] ap- 
proach for analysis boson eigenfrequencies. We also note 
that the present study can be extended to take into ac- 
count finite temperature effect [s^l on the excitations. 
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